Performance of LCD iterative method in the nite element and nite di erence solution of convection–di usion equations
نویسندگان
چکیده
In this work, we evaluate the performance of the left conjugate direction method (LCD) for the solution of non-symmetric systems of linear equations arising from nite element and nite di erence discretizations of the convection–di usion equation. We extend the LCD algorithm proposed by Dai and Yuan (Int. J. Numer. Meth. Engng 2004; 60:1383–1399) to accommodate restarts. Our discussion considers comparison studies between the computational e ciency of the GMRES and LCD methods and some issues related to the choice of the forcing term in the inexact Newton method. Copyright ? 2005 John Wiley & Sons, Ltd.
منابع مشابه
Finite Volume Methods for Convection Diffusion Problems
Introduction In this paper we consider cell centered nite di erence approximations for second order convection di usion equations of divergence type Our goal is to construct nite di erence methods of second order of approximation that satisfy the discrete maximum principle The error estimates are in the discrete Sobolev spaces associated with the considered boundary value problem Approximation ...
متن کاملLeast squares streamline di usion nite element approximations to singularly perturbed convection di usion problems
In this paper we introduce and study a least squares nite element approximation for singu larly perturbed convection di usion equations of second order By introducing the ux di usive plus convective as a new unknow the problem is written in a mixed form as a rst order system Further the ux is augmented by adding the lower order terms with a small parameter The new rst order system is approximat...
متن کاملHybrid Laplace transform ®nite element method for solving the convection±dispersion problem
It can be very time consuming to use the conventional numerical methods, such as the ®nite element method, to solve convection± dispersion equations, especially for solutions of large-scale, long-term solute transport in porous media. In addition, the conventional methods are subject to arti®cial diusion and oscillation when used to solve convection-dominant solute transport problems. In this ...
متن کاملA Robust Algebraic Preconditioner for Nite Diierence Approximations of Convection{diiusion Equations
Stable nite di erence approximations of convection-di usion equations lead to large sparse linear systems of equations whose coe cient matrix is an M{matrix, which is highly non symmetric when the convection dominates. For an e cient iterative solution of such systems, it is proposed to consider in the non symmetric case an algebraic multilevel preconditioning method formerly proposed for pure ...
متن کاملPreconditioning a mixed discontinuous finite element method for radiation diffusion
We propose a multilevel preconditioning strategy for the iterative solution of large sparse linear systems arising from a nite element discretization of the radiation di usion equations. In particular, these equations are solved using a mixed nite element scheme in order to make the discretization discontinuous, which is imposed by the application in which the di usion equation will be embedded...
متن کامل